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Abstract 

A personal summary of the discussions which took place at the 
informal meeting in Amirim, Israel from June 1-4 2000, concerning 
the dipole picture of small-x physics is presented. The broad aim 
of the meeting was to address the question "Has HERA reached a 
new QCD regime (at small x) ?". The new regime in question is the 
high-density, but weak-coupling, limit of perturbative QCD. 



1 Introduction and background 



From 1-4 June 2000 approximately twenty HEP experimentalists and the- 
orists met at Amirim, a small vegan village in Northern Israel, to discuss 
the question "Has HERA reached a new QCD regime (at small x) ?". The 
meeting was organized by Tel Aviv University and was somewhat unusual 
and innovative in style in that discussions took place "al fresco" , at a white 
board in a small garden, and there were no formal presentations. This paper 
aims to summarize the discussions and is written with the assistance of those 
who attended the meeting^ 

In recent years improvements to the HI and ZEUS detectors, have enabled 
the measurement of the Deep Inelastic Scattering (DIS) structure function, 
F2{x,Q'^) to be pushed to new and interesting kinematical regions in the 
photon's virtuality, Q^, and scaling variable x ~ Q'^/W'^ {W is the centre 
of mass energy of the photon-proton interaction). In particular, the new 
high-precision very low <^ 1.0 GeV^ data ^ |^ allows the important 
transition from DIS to photoproduction to be studied, at low values of x, for 
the first time. It is known that F2 must be proportional to as — *• (to 
satisfy conservation of electromagnetic current and to achieve a finite pho- 
toproduction cross section) and indeed this behaviour appears to have been 
observed by ZEUS at the lowest values (see Fig.(|l]) and |[l|). This implies 
a departure from the logarithmic scaling violations predicted by standard 
perturbative QCD evolution analyses (leading- twist DGLAP[^) which are 
anyway not expected to be applicable for such low scales (the input scale Qq 
is currently chosen to be around 1.0 GeV^ and all data with > Qq are 
analyzed, implying that the higher-twist corrections to the standard picture 
are not important above Qq). 

To illustrate this transition region, the slope dF2{x,Q'^)/d\n{Q'^) as a 
function x is often plotted in the low region, to search for deviations 
from the logarithmic scaling violations predicted by perturbative QCD. This 
derived quantity is sensitive to the gluon density at small x. Figs.(0J^) show 
the latest plots of the slope derived from the HI and ZEUS data at fixed 
values of (see e.g. ^ for more details). In the measured range, the 

slope is unambiguously seen to increase monotonically as x decreases even 
for as low as 0.75 GeV^. An earlier, cruder version of this plot appeared 
to show a turnover in x |^ and created a great deal of interest and discussion. 
It now appears from Figs. (H, ^ that this was a result of the limited range in 

accessible at a particular x, and a changing range at different x values, 
in the older data. 

*See Acknowledgments. 
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The lower plot of Fig.(^) [0] shows the slope for fixed W, rather than Q^. 
A genuine turnover in x is revealed when the data is presented in this way. 
The upper plot shows the slope plotted as a function of Q^, again revealing 
a turnover for fixed W. It is interesting to note that the position of the 
maximum appears to shift to higher as W increases. 

The small-x limit, which is also probed at HERA, corresponds to the 
high energy, or Regge, limit of DIS {W'^ ^ Q'^^nip, ...). The convergence of 
QCD perturbation theory in as(<5^) is known to break down eventually here 
due to logarithms in 1/x (see e.g. and references therein). Although the 
standard global DGLAP QCD fits produce an adequate description of the 
inclusive F2{x,Q'^) data, over a wide range of {x,Q'^), certain characteristic 
features of the resulting fitted parton densities at low and low x suggest 
that this approach may require significant corrections. For example, the 
DGLAP analysis translates the experimentally-observed steep rise of F2 at 
small X into a numerically large gluon density, which rises steeply with 1/x 
(cf. Figs.(H, II)). Clearly, if parton densities get too large the partons must 
start to interact with one another in order to tame this growth (see e.g. [0), 
leading to a breakdown of the standard DGLAP approach. This corresponds 
to the high density but weak coupling limit of perturbative QCD. 

In addition, a startling experimental observation at HERA concerning 
DIS in the small-x region is that an approximately constant fraction (of the 
order of 10%) of events for x < 10^^, are diffractive (i.e. the proton remains 
intact rather than scattering inelastically). This is rather difficult to explain 
from within the standard parton model picture of a 'typical' DIS event: in a 
'fast proton' frame, a highly virtual photon strikes a charged coloured quark, 
carrying collinear momentum fraction x of the proton. This violent collision 
generates strong colour forces that rip the proton bound state apart. In this 
scenario, large rapidity gaps in the final state, which characterize diffractive 
events experimentally, are expected to be exponentially suppressed. 

In the past few years the theoretical description of small-x DIS has shifted 
somewhat from this 'infinite momentum frame' to the frame in which the pro- 
ton is at rest. This is particularly useful for considering diffractive scattering 
(and especially the subclass of diffraction in which vector mesons and real 
photons are produced exclusively). In this frame, in the high energy limit, 
there is an appealing factorization of small-a; processes which in some sense 
is purely kinematic: the virtual photon fluctuates into a hadronic system at 
large distances upstream of the proton target, this hadronic state takes a 
long time to form, interacts briefly with the target and subsequently forms 
the hadronic flnal state over a long period. From this frame it is much easier 
to believe that if the interaction is "gentle" the integrity of the proton bound 
state will not be compromised. 
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The simplest fluctuation, which dominates for very small transverse size 
systems is the quark anti-quark state which forms a colour-triplet dipole. 
Because of this, the generic name of "the dipole picture" or "the dipole 
approach" has been adopted for proton rest frame descriptions. For larger 
transverse-size systems more complicated fluctuations, containing more than 
two partons, must also be considered. However, even for such complicated 
systems the transverse diameter of the fluctuation impinging on the target 
remains a useful variable to classify interaction strengths of the fluctuations 
of the virtual photon. For linguistic convenience, the transverse diameter and 
the interaction strength have become synonymous with the misnomers dipole 
size and dipole cross section respectively, regardless of the fact that more than 
two partons are typically involved for larger transverse-size fluctuations. 

Within this framework, inclusive cross sections, for example ai, may be 
written in the following, generic, factorized way []TT], |T^, |T^: 



1 
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aL{x,Q ) = J dz J d d±a{dj_, ...)\ilj^^L{z,d±)\ . (1) 

In this simple case, the light-cone wavefunction of the longitudinally polarized 
virtual photon, in terms of the qq is given by QED: 

6 ""^ 

\M^,d^)\' = J]e2QV(l - zfK',{ed^) , (2) 

^ 9=1 

where = Q'^{z{l — z))+m'^ and z is the fraction of the photon's momentum 
carried by the quark (of mass rrig). The variable d± is Fourier conjugate to the 
relative transverse momentum of the gg-pair and in general corresponds to the 
transverse size (diameter) of the scattering system. An appealing feature of 
this representation is that the short and long distance contributions to a given 
process are manifest. The power of the approach lies in the possibility that 
the interaction strength depends only on the configuration of the interacting 
system and not on how it was produced. For example the strength depends 
on the system's transverse size and not on external variables such as Q^. 
This implies a universal interaction cross section, a, which may be used in a 
wide variety of small- a; inclusive, exclusive and diffractive processes. 

There are several models available for this interaction cross section a 



and we compared and contrasted four of them [|T^, |T5|, [T^, in detail at the 
meeting. A summary of these discussions is given in the Appendix. One issue 
is precisely how this cross section depends on the energy of the interaction, i.e. 
whether it depends on or x (bearing in mind that it must reproduce the 
observed approximate scaling properties of the structure function F2{x, Q"^)). 
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In principle, a can also depend on z, although all models constructed so far 
neglect this dependence. There are many approaches influenced by similar 
ideas covering various aspects of small-x physics. In this short note we make 
no attempt to be complete in covering all references and instead refer the 



reader to the references within 14, 15, 16, 171 and to recent review articles 



on diffraction [OB, m 20 



Broadly speaking, each of the models discussed share some common fea- 
tures designed to encapsulate the important physics issues. For small ci^ the 
dipole cross section is ascribed a steep rise with energy (from the point of 
view of QCD this corresponds to the steeply rising gluon density). A much 
more gentle (or even flat) behaviour in energy is ascribed for larger non- 
perturbative sizes (in agreement with observations of hadronic cross-section 
data). There are two important consequences of this. Firstly, each model 
contains a mixture of perturbative and non-perturbative aspects. The inter- 
play between them may be systematically studied, within a particular model, 
by examining the (i^-integral for the process of interest. Generally, the hard 
perturbative part becomes more important as the energy increases (since it 
grows faster than the soft piece). Secondly, the cross sections associated 
with small dipoles reach hadronic sizes (tens of millibarns) in the kinematic 
region probed at HERA (because the gluon density is numerically large at 
small x). The latter feature implies that some kind of "taming" or satura- 
tion physics is necessary at high energies in order to preserve the natural 
geometrical criterion that the dipole cross section increases monotonically 
with transverse size. The introduction of such taming effects implies a new 
high-density regime of hard QCD interactions. 

Our primary goal is to identify clear signs for taming effects in experimen- 
tal measurements. Of course, in order to claim that we have found something 
new, we need to understand more precisely how the dipole picture relates to 
the more rigorous DGLAP formalism in the kinematic region where they are 
supposed to agree (moderate and small, but not too small, x). In prac- 
tice, since the analysis of hadronic interactions is always complicated by our 
poor understanding of confinement, it will probably be necessary to build up 
a body of experimental evidence, combining data from several experimental 
channels, in order to force the unambiguous conclusion that we have entered 
a new regime. 

The connection of the dipole approach with standard DGLAP is con- 
sidered in Section 2. The issue of 'higher twist' corrections to DGLAP is 
discussed briefly in Section 3. Section 4 covers the description of inclusive 
diffraction and exclusive diffractive processes (such as vector meson produc- 
tion) within the dipole framework. The connection to the DGLAP evolution 
of diffractive parton densities is also briefly discussed. Finally, in Section 5 
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we conclude with a 'wish hst' for better measurements of particular exper- 
imental processes that will pin down the details of the dipole cross section 
and may reveal saturation phenomena, thus achieving a principle aim of the 
meeting, i.e. to help prioritize the experimental effort in this area. 

We note in passing that although we were mainly concerned with high en- 
ergy interactions of photons and protons, the concept of a universal high en- 
ergy interaction cross section, embodied in a, is in principle more widely ap- 
plicable. It may be used for example for DIS using nuclear targets, in photon- 
photon collisions and even for certain hard processes in hadron-hadron inter- 
actions. 

2 Comparison of the dipole picture with the 
standard DGLAP analysis 

Clearly the dipole approach to small-x processes should be complementary to 
the DGLAP analysis (for moderate and not too small x). In this section 
we briefly recap the standard approach and then emphasize similarities and 
differences between the two pictures. 

The standard perturbative QCD approach to calculating the structure 
functions of DIS invokes collinear factorization and involves a convolution of 
universal long-distance parton (quark and gluon) density functions ( "pdfs" , 
into which all infrared divergences are factorized) specified at the factoriza- 
tion scale, /ij, with short-distance, hard coefficient functions for the fusion of 
these partons with the virtual photon. The coefficient functions are calcula- 
ble in perturbative QCD as a power series expansion in as- This power series 
generates logarithms in to all orders in perturbation theory (usually 

one chooses /Xj = to suppress these potentially large logs in the coefficient 
functions and "move them" to the pdfs). The independence of physical cross 
sections on the choice of this factorization scale then leads to the DGLAP 
equations which embody the physics of the renormalization group and 
describe how the pdfs change as their scale changes. The pdfs, which are 
solutions to the DGLAP equations, are then specified to a given order in 
logarithms (LO, NLO, NNLOQ,--.), with the accuracy increasing as more and 
more terms are included in the perturbative expansion. The parton densities 
at the input scale are not calculable in QCD, but factorization implies that 
they are universal functions. In practice, the various groups responsible for 
producing pdfs (CTEQ, MRS, GRV) constrain their analytic forms in x at 

^^For recent progress on NNLO calculations for parton densities and structure functions 
see e.g 0. 
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the input scale, Qq, using sum rules, symmetries, etc. The non-perturbative 
parameters used to specify the input distributions are then determined by 
performing a minimum fit to all available relevant data, evolving to the 
appropriate scale using leading-twist DGLAP evolution. The value to choose 
for the input scale, Qq, is not predicted by the theory. Recently scales as low 
as Qo = 1-0 GeV^, or even lower in the case of GRV [^], have been used. 



The factorization is termed collinear because all transverse momenta (de- 
fined relative to the photon-proton axis) in the pdfs is forced to be less than 
li'j. As far as the photon is concerned it strikes a parton which has negligible 
transverse momentum relative to the proton's direction, so all input partons 
are moving collinear with the proton. 

From the theoretical point of view this collinear factorization, which sep- 
arates short and long distance contributions, is made possible by the use of 
the operator product expansion (OPE). For very large Q^, a subclass of all 
possible Feynman graphs are dominant in which, for a particular choice of 
gauge, the hard coefficient function is connected to the proton by only two 
parton lines (the so-called "handbag graphs" and their perturbative QCD 
corrections). Clearly, for sub-asymptotic this picture has corrections 
which the OPE predicts contribute at relative order (9(1/Q^) and beyond 
{0{1/Q'^)"',n = 2,3,...). These are commonly called "higher-twist correc- 
tions", although it should be remembered that within the OPE twist has 
a rigorous definition^. For example, a subclass of the most important sub- 
asymptotic corrections (twist four) corresponds to the set of Feynman graphs 
in which the hard coefficient function is connected via four parton lines to 
the proton. 

For very small x the standard implementation of DGLAP must break 
down due to the development of logarithms in energy or 1/x at each order 
in perturbation theory, which eventually must spoil the convergence of the 
perturbative series because the product ^^(Q^) ln(l/a;) becomes C(l). The 
leading logarithms in 1/x are re-summed by the BFKL equation and 
recently next-to-leading logarithmic corrections to the leading-log result were 
completed |^ and were found to constitute a very large correction to the 
leading-log result. Subsequently, considerable progress has been made in 
understanding the reasons for the large size of these corrections (for a recent 
overview and references see [p5[1). 

A long standing objection to the use of BFKL calculations was again 
raised by L. Frankfurt in the discussion session, i.e. that the HERA kinemat- 
ics provide a relatively restricted rapidity range for perturbative real gluon 

^FoT a given operator the twist is the difference between its mass dimension and spin, 
and for a given product of operators it is the sum of the twists of the individual operators. 
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emission in quasi-multi-Regge kinematics. Realistically, HERA kinematics 
allow for typically two or three real gluon emissions on average, so that sum- 
ming an infinite tower of gluon ladders may be unrealistic. Moreover, the 
application of these re-summed calculations to relevant physical processes 
are plagued by many uncertainties, including in many practical situations 
an uncontrollable contamination of non-perturbative physics. The final so- 
lution to this extremely difficult problem is far from clear. This fact in itself 
lends support to, and makes timely, phenomenological frameworks such as 
the dipole picture which embody a certain amount of modelling concerning 
the mixing of perturbative and non-perturbative physics. The interplay of 
new experimental results and the refinement of these models then becomes 
of paramount importance in attempting to pin down the relevant physics. 

For small x it is necessary to generalize coUinear factorization to what 
is known as "fc^-factorization" [^] and to discuss unintegrated distributions 
rather than conventional pdfs (the subject of the BFKL equation is the un- 
integrated gluon density). At leading- twist level, the former are related to 
the latter via an integration over transverse momentum, e.g. for gluons 



x9{x,Q')= (3) 
Jo '± 

The connection between the dipole picture and standard DGLAP (and es- 
pecially /cx-factorization) is far from obvious. Clearly both pictures involve 
a photon connected to a quark loop. The transverse momentum and en- 
ergy fractions used to specify the ipi^, k±)s in the dipole approach are closely 
related to the momentum fiowing around this loop. It would be extremely 
interesting and useful from a theoretical point of view to understand collinear 
factorization in d^-space more precisely. 

One recent model claims to include as much information as possible 



from DGLAP in that it is explicitly built using gluon densities from DGLAP 
fits (implicitly incorporating DGLAP evolution) and the well-known leading 
In formula for the dipole cross section: 

cTpQCDix, d\) = —^as{Q'^)xg{x, Q^) , (4) 

which has been widely used in the analysis of various diffractive channels (for 
a review see Ijl^). Indeed, two-gluon exchange has been explicitly proven to 



be the dominant contribution in certain hard exclusive processes p7| , p8 
The eikonal model of [T^ also reduces to precisely this form for the limit of 
very small d± (see Appendix). 
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Unfortunately, an additional assumption concerning the connection be- 
tween transverse sizes and four-momentum scales is unavoidable. This leads 
to a specific model-dependent assumption about the connection between the 
dipole picture and standard DGLAP. Having made this connection (using 

= 10 /d\) the model of 0] uses eq.@ over as wide a range of transverse 
sizes as possible, i.e. for all dj_ < dj_^Qo (where (i_L,Qo = VlO/Qo corresponds 
to the input scale Qo). The model of flSl makes the choice d±^Qg = 2/Qq 
instead. Of course, both choices are equivalent to leading-log accuracy. 

The nature of evolution in the dipole picture is particularly interesting. 
At a particular Q"^, the integrand of a given cross section, such as aL of 
eq.(|lD, has a profile in d±. In this sense, the d±-pTofi\e for F2 at = 
Qq "corresponds" to a model for the input quark density. The light-cone 
wavefunction of the photon clearly changes as increases and this effect 
"evolves" the profile to smaller d±. This also means that the gluon density 
in eq.(|D is predominantly being sampled at typically larger scales. This 
shift in the profile in d± with is very closely related to standard DGLAP 
evolution but unfortunately, a more precise connection can not be made at 
present. 

The other two models |16|, were "designed" to describe the transition 
from the photoproduction region to DIS (from small to moderate Q"^). Hence, 
the known QCD behaviour at small x (steeply rising gluon densities) was 
included in a more phenomenological way by including a piece that rises 
steeply with energy at small dj_ (see also [^). 

Overall at the meeting there was a broad consensus that eq.(|) imposes 
the correct QCD behaviour, to leading-log accuracy, for small dipole sizes. 



3 Higher twists 

Unfortunately, from a theoretical point of view, very little can be said about 
the numerical importance of higher twist effects, at a given Q^, with a high 
degree of rigour. 

One of the most basic and perhaps one of the more important higher 
twist corrections involves graphs where four gluons couple directly from the 
proton to the quark box. An analysis |30| of such graphs revealed an intrigu- 



ing partial cancellation between the correction to the longitudinal structure 
function, Fl, and the transverse one. Ft, when they are summed in F2 (see 



also for an analysis of related graphs). This cancellation has been taken 
by some supporters of the standard approach as an explanation as to why the 
leading-twist DGLAP analysis is so successful. In addition, this cancellation 



was also seen explicitly within the context of a specific model |T^ for a. If 
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the cancellation is genuine in F2 it provides yet another strong motivation 
to make an explicit experimental measurement of in the region of moder- 
ate (for example by lowering either the proton or lepton beam energy at 
HERA) to look for deviations from the leading-twist predictions. Of course, 
in order to be discriminatory, this would require a sufficiently long running 
period to make a high precision measurement of Fl. 

During the discussion sessions the point was made by L. Frankfurt that 
conclusions based on specific subclasses of twist-four graphs may be too naive. 
He speculated that perhaps once twist-four corrections become numerically 
important all other higher twists may as well (i.e. the twist expansion could 
break down). This scenario would imply that the system undergoes some 
kind of phase transition to a completely new regime. 

4 Diffraction and heavy flavours 
4.1 Exclusive processes as probes of a 

In hard exclusive processes, such as diffractive vector meson production and 
deeply virtual Compton scattering, a enters on the amplitude level, rather 
than on the cross-section level. For example, the differential cross section for 
vector meson production is given by: 



where ipv{z,k±) is the light-cone wavefunction of the produced vector me- 
son. Thus, these exclusive processes are sensitive to integrals that weight a 
differently in d± than inclusive structure functions. As an example, Fig.(^) 
shows the integrand, or profile, in d± of the amplitude for photoproduction 
of J/ip at three different values of W in the HERA range. A hybrid light- 
wavefunction for J/ ip from was used and the model of |[T^ for a. Moving 
to larger Q'^ would cause the peak and profile in d± to shift to small distances 
as the process becomes more perturbative. Varying and and looking 
at different states (7, p, (p, J/ip, T) allows a detailed scanning of a over a wide 
range in d± and energy. Therefore, in principle exclusive processes provide 
complementary information to inclusive structure functions that may be used 
to pin down the detailed functional form of the dipole cross section. 

Unfortunately, these processes are sensitive to several other uncertainties, 
including the effects of skewedness of the amplitude, model dependence asso- 
ciated with ^y(^, k±), poorly known t-slope, etc, which make uncovering the 



da 
It 



t=o oc \A{s,t = 0)W 



(5) 



A 
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details of a a difficult job in practice. To obtain the maximum information, 
a simultaneous analysis of all available data within the same framework is 
required. 



4.2 Inclusive diffraction 

One of the appealing features of the dipole picture is its ability to describe 
diffractive and inclusive scattering within the same framework. The qq dipole 
produced by a longitudinally-polarized photon has the following contribution 
to inclusive diffraction: 



da 



dt * ° 167r 



d'd^dz\^lj^^L{z,duQ )m ■ (7) 



Note that a appears squared. Roughly speaking, one may think of this as 
the square of the amplitude in eq.®, but instead of an explicit vector meson 
wavefunction appearing completeness in the sum over all possible produced 
states is used. 



The proof by Collins of factorization in hard diffractive scattering ||33 
placed the DGLAP analyses of diffraction in DIS on a ffimer theoretical 
footing. The fundamental objects which evolve with are known as diffrac- 
tive parton densities. Relative to standard pdfs, they are defined with the 
additional constraint that the proton remains intact in the final state. As 
such they depend on two extra variables which parameterize the scatter- 
ing proton (a;jp, the longitudinal momentum fraction lost by the incoming 
proton, and t, the square of the four-momentum transfered to it, are usu- 
ally chosen). With the additional assumption of Regge or Ingleman-Schlein 



factorization, a DGLAP analysis [Q, ^ leads to the conclusion that the 
diffractive gluon density is much larger than the diffractive quark density for 
values of/? = x/xjp >0.1 (the ratio of gluons to quarks is several times larger 
in diffraction than in the inclusive case). 

In diffraction, the connection between the dipole picture and DGLAP- 
driven diffractive parton densities is more obscure. Qualitative statements, 
similar to those made in the inclusive case, concerning the evolution of the d±- 
proffie as changes may also be made here. A further essential complication 
is that one must also include higher order Fock states (including qqg) in the 
virtual photon, especially for the case of large diffractive masse^. The ratio 
of diffractive to inclusive total cross sections is particularly interesting since 



^In contrast, inclusive structure functions are unique in that it is possible to factorize 
all gluon emission into the definition of the dipole cross section. For hard vector meson 
production the qqg contributions are suppressed by the light-cone wavefunction of the 
vector meson. 
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it also helps to specify the dominant d±. The fact that this ratio appears to 
be fairly independent of energy supports the widely-held view that inclusive 
diffractive processes are much more sensitive to larger transverse sizes. 



4.3 Diffractive heavy flavour and jet production 

Diffraction can also occur for processes involving a colour-octet pseudo-dipole 
in the initial state. These result from asymmetric qqg fluctuations in which 
the quark and anti-quark remain close together in transverse space. This 
may be achieved by demanding heavy flavours or a large relative p± between 
the quark and anti-quark, with the gluon carrying a small fraction of the 
available photon energy. Such fluctuations are closely related to diffractive 



gluon densities (see [|T8[ for more details) and have become known as the 
"gluon channel" for diffraction. These configurations are known to have 
a large colour factor (9/4) relative to the colour-triplet dipole (the "quark 
channel"). This large factor may make the overall a for hard gluon-channel 
diffractive processes much closer to the saturation limit. As such, it should 
provide a clear signal for new physical effects in the saturation regime (see 
p6| for a longer discussion of these issues). The difficult part is finding a 
physical process in which one is sensitive to small enough colour-octet dipoles 
to observe taming effects in a perturbative environment. 

If the large colour factor is translated directly to a it leads to some strik- 
ing predictions for diffractive channels that are driven by the exchange of 
gluons and dominated by colour-octet fluctuations of the photon. One excel- 
lent channel to consider is diffractive heavy flavour production. The above 
observation, concerning the large ratio of diffractive gluon density to diffrac- 
tive quark density relative to its inclusive counterpart, naively suggests the 
striking prediction that about 30 — 40% of charm quark production should 
be diffractive. Realistic kinematic constraints, including a typical effective 



scale of about 10 GeV^, implemented in the RAPGAP Monte-Carlo P? 



suggest that in practice one may expect a ratio closer to 20%. This effect 
is not yet seen in the currently measured region, which is concentrated at 
rather low values of /3 ~ Q^/ {Q^ + M^) because of the necessity to produce 
the heavy charm quark pair (according to ZEUS the diffractive charm 



appears to be about 6% of the total). With the expected higher luminosity 
it will be possible to collect more data at higher values of and hence (3. 
An observation of a sharp increase in the diffractive charm ratio would be a 
strong indication of the correctness of this picture. 

Another gluon-driven channel is the diffractive production of jets (both 
two and three jet events are interesting). One experimental problem here 
is that the requirement to have a well-defined jet naturally forces one to 
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consider a large Mx system and hence one is pushed again into the low (3 
region (high a;jp at fixed x). Again, the higher luminosity will allow more 
data to be collected at higher (3. 

5 Conclusions and a wish list of desirable mea- 
surements 

The HI and ZEUS experiments have each collected over 100 pb"^ of data 
so far. With the increase of luminosity after the HERA upgrade comes 
the opportunity to improve considerably the precision of a large number of 
measurements of diffractive processes. At this meeting, we tried to prioritize 
which of these processes are the most important to measure well in order to 
determine whether a new saturation regime really has been reached at small 

X. 

From the point of view of the dipole picture we need better measurements 
to pin down the energy and transverse size dependence of the universal func- 
tion a. As discussed above, this may be realistically achieved only by a 
combined analysis of several exclusive and inclusive processes. 

The precision measurements of F2 are now excellent and very high pre- 
cision data exist over a very wide kinematic range. Indeed, the main uncer- 
tainties in F2 over a wide range in are now given by systematic rather 
than statistical errors. It would therefore be extremely useful to have the 
complimentary information that a dedicated high-precision measurement of 
Fi{x,Q'^) in the same kinematic range would provide. This would help to 
constrain the dependence of a on d\_ and with sufficient precision could act as 
a test of the DGLAP parton densities. A serious deviation from the DGLAP 
predictions (determined largely from the measurement of F2) may confirm 
the significance of higher twists in F^ (see section 3). 

Improved statistical precision on the derived quantity dF2/ d\n.Q^ in the 
region of relatively low would be a good discriminator between approaches 
and a combined quantitative analysis of this information with the description 
of J/ip production within the same framework was advocated by U. Maor. 
He claimed that none of the currently available pdfs were able to fit both 
quantities simultaneously, without including taming effects []. Of particular 
interest is further constraining the energy dependence of hard (either high 

or large mass My) vector meson production in the highest possible bins 
in W"^. An observation of a flattening of the steep rise with energy would 
constitute a clear and relatively model-independent signal for taming in these 

^Since the meeting a preprint ||39| has appeared. 
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processes. 

Improved measurement of diffract ive charm and diffract ive dijets, espe- 
cially at higher values of would be particularly welcome since they are 
sensitive to gluon induced processes in which saturation effects may be seen 
early. 

There was a consensus of opinion that it was vital that the t— dependence 
of various diffractive processes should be measured, including, if possible, how 
it changes with energy (shrinkage). The t-dependence is usually parameter- 
ized in exponential form, e'-^^*-', using the t-slope parameter, B£,{yV'^\ The 
rate of decrease of this diffractive peak with increasing t and W'^ for a given 
diffractive process, embodied in gives a strong indication of the mixture 
of soft and hard contributions. It was agreed that the measurement of Bd is 
of major importance. 

From an experimental point of view, we are strongly encouraged by the 
proposal of HI to insert a dedicated detector in the very forward region to 



detect scattered protons This detector will allow a clean tagging of 

diffractive interactions, taking full benefit of the high luminosity in the new 
HERA running regime. This is particularly important for measurements of 
the diffractive structure function (where the major uncertainty is associated 
with the presence of proton dissociative background), and of jet and of charm 
diffractive production. Most relevant for the present discussion will be the 
possibility of measuring the slope parameter, By). 

In conclusion, the first hints of a new saturation regime may already 
have been seen at HERA, in a rather indirect way, using the conventional 
framework of DGLAP. The phenomenological dipole approach is capable of 
producing clear predictions for experimental signatures that would reveal this 
new regime and the relevant calculations are already underway. There are 
excellent prospects for improved measurements of many small-x processes, 
which are capable of confirming this new regime directly. In practice, it 
appears to be necessary to perform a combined quantitative analysis of cur- 
rent and future data for several different channels to build up the necessary 
evidence. 



Appendix: Features of different models for (j 



The aim of this section is to compare and contrast several models [14, 15, 



16| , p!7| for the dipole cross section and how they depend on and energy. 
This section is included mainly for information and no attempt is made to 
be critical of physics used in building particular models. A more critical 
discussion from the point of view of the author and his collaborators may 
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be found in Section 4 of While much of the discussion in Amirim, 

perhaps inevitably, centred on the contrast between the different approaches, 
a tentative consensus of the discussion sessions was that the features which 
the various models share should eventually prove to be much more significant 
than the details in which they differ. The most important of these shared 
features was the implementation of some kind of taming of the rapid growth 
of the short distance part of a as the energy increased. 

First of all we consider the simple saturation model for a suggested by 
Golec-Biernat and Wiisthoff |[T6[| : 



cr(x, dl) = ao(l - eM-dlQl/A{x/x^Y]). (8) 

A three parameter fit to the HERA data on F2 with x < 0.01, excluding 
charm and assuming Qq = 1.0 GeV, produced the following values do = 
23.03 mb, xq = 0.0003, A = 0.288 and a reasonable x^- This model also 
did an adequate job in describing the diffractive data (see second reference 
of |T^) without any further tuning of parameters (but with some additional 
assumptions, for example that exactly the same a also couples, in a particular 
way, to the qqg component which is required to describe inclusive diffraction). 

Certain features are apparent in eq.(^. First of all, a is assumed to de- 
pend only on the transverse size, d±, and upon the Bjorken scaling variable 
X. For small transverse size, -C 4(a;/a;o)'*'/Qo) ^ designed to be propor- 
tional to the area of the scattering system (a feature expected on geometrical 
grounds) and to grow steeply with 1/x at small x: a{smal\ d±) oc d^ x~^'^. In 
contrast, at large d± flat behaviour in d± and x is assumed: (T(large d±) ao. 
As X decreases this fiat region extends to smaller and smaller d±. 

For certain exclusive electroproduction processes such as Vector Meson 
Production and Deeply Virtual Compton Scattering |^ the interaction 



of a small qq has been rigorously proved to proceed predominately through 
the exchange of two gluons in the t-channel. Thus perturbative QCD relates 
a, to leading \og{Q^) accuracy, to the leading-log gluon density in the proton 
(for a detailed derivation see ^T|): 



TT^d^ 

^pQCd{x, d\) = —^asiQ'^)xg{x, Q"^) . (9) 

To leading-log accuracy one may make the choice = Q^. 

A detailed model for a based on this form is given in fl^ for all transverse 
sizes. The detailed form relies on the continuity and monotonic increase of 
(T as a function of dj_ (at fixed x). For large, non-perturbative, transverse 
distances a is matched on to a typical hadronic cross section (pion-proton is 
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used at d±_ = d±^Tr = 0.65 fm). This is ascribed a slow Donnachie-Landshoff 
energy growth typical of hadronic cross sections fl 

The model in attempted to go beyond leading-log accuracy in a by 
incorporating a (i^-dependence for the four- momentum scale = 10 /d\. 
In addition, momentum conservation is implemented for the fusion of the 
photon and the gluon. This leads to the gluon being sampled at an effective 
momentum fraction x' which depends explicitly on the mass of the produced 
gg-pair (which may be translated into a d± dependence in the simplest case). 
With these modifications in place the form given in eq.(^ is used within its 
region of applicability, i.e. for d^ < 10/ Ql, providing x is not too small. Note 
that the generic geometric d^ behaviour becomes modified by the behaviour 
of OgXg on which translates into a complicated effective behaviour on d± 
(which will depend on the value of x). 

For very small x, it was found that using CTEQ4L leading-log partons 
the numerical value of the dipole cross section at small d± became large (of 
hadronic interaction strength) due to the numerical size and steep growth 
of the gluons at small x. This observation threatens to spoil the monotonic 
increase of a with dj_ expected on geometrical grounds and was taken to 
be a signal for the necessity of taming corrections to the exchange process. 



These were implemented by hand in [|T4| at fixed x at the critical point, 
d±,crit{x), at which a reached 50% of its maximum value (as determined by 
the pion-proton cross section). The perturbative QCD expression for a in 
eq. (^ then has a progressively reducing region of applicability as x decreases: 
d± < dj_^crit{x) (cf. the "critical line" of |]16[ and the packing factor of [|15]). 



As an example of a comparison between models, Fig.(|g) shows a plot of 
the ratio of these two models for a. The two forms have inherently different 
behaviour in energy and transverse size which is especially apparent at small 
and large d± in HERA kinematics. For intermediate d± the numerical values 
for a are similar. 

Forshaw, Kerley and Shaw [|T^ proposed a general ansatz for a, which 
they stress is modelled as a function of d^ and (rather than x), with a 
soft Pomeron and a hard Pomeron piece: 

d^)=a ^ 5pf^J ^ (^^^')'" + dlPl{dM-''l'-\dlwY'^ , (10) 
1 + P^{di_) 

where Ps{di_) and Ph{di/) are polynomials in d^. They successfully fit this 
general form to the inclusive data. For very large dipoles (of the order of the 

"in order to preserve the uniqueness of ct as a function of x and dj^, it was necessary 
to implement this energy growth via a multiphcative factor (x/xo)'^ '^^ rather than the 
usual (s/so)'^ '^^ behaviour. Since xq and sq are in principle unknown, this is an acceptable 
compromise over a reasonable range in energy. 



16 



size of a light meson) one may think that the dipole cross section should be 
a function of rather than x ~ Q'^/W^, since it should just depend on the 
energy of the collision (here there is no hard scale with which to specify an 
'x'). However, in order to achieve the observed approximate Bjorken scaling 
of F2, it was necessary to model the dipole cross section in eq.([lO|) as a func- 
tion of the dimensionless d\W^, rather than alone. This model also gave 
a reasonable description of the diffractive data (see second reference of [|T7| ) 
without further tuning of parameters, but with some additional assumptions 
on the qqg component. 

Total cross sections are related to the forward {t = 0) part of the elastic 
scattering amplitude via the optical theorem. However, the amplitude itself is 
of course a function of s and t = n^n^ ^ — ra^ (in high energy kinematics the 
momentum transfer is mainly transverse). In general, one may also consider 
the Fourier transform with respect to this momentum transfer: 



A{s,t) = J d'b^e-'''^'^A{s,b^) . (11) 

For t = the amplitude can be represented as an integral over all impact 
parameters b±, 



A{s,t = 0) = j d'^b^A{s,b^). (12) 



More generally, the cross section for a given process will also depend on 
the impact parameter of the collision, a quantity that depends on both the 
scattering state and the target and should not be confused with the transverse 
diameter of the dipole, d±. 

Several years ago Gotsman, Levin and Maor [|15| developed a form of the 
eikonal model to be applied to DIS processes 

a{x,d^) = j d%^e-'''^'^^a{x,d^,b^), (13) 

a{x,d^,b^) = 2(l-e-^/2). (14) 
The following convenient factorized form was assumed for the eikonal 

71-2 4 

n = S{bi_) — d\ as xg{x,-^) , (15) 

with a Gaussian form ascribed to the normalized non-perturbative two-gluon 
form factor of the target: 



"-AT 



*We follow closely here the presentation in Section 2.2 of the second reference of nl 
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where R% is related to the proton radius. It is this specific Gaussian form 
in b± which leads to the popular exponential parameterization of the t- 
dependence, e^^*. In the limit of small Q, the exponent of eq.(|Tip may 
be expanded and the familiar form of eq.(^ is recovered for a. As stated 
earlier, this corresponds to the exchange of a pair of gluons in the t-channel, 
treated in the leading log((5^) approximation and constitutes the basic build- 
ing block in this version of the eikonal model. In determining whether Q is 
indeed small one must examine the magnitude of the packing factor of par- 
tons in the cascade, defined as 

_ ^pQCD{x,d±) 

^ " 7rR% ■ ^^^^ 



The line k = 1 is analogous to d±^crit{x) of |jT^ and the critical line of |T^. For 



large fi, a standard Regge analysis is adopted. Thus the energy dependence 
of this model is similar, but not identical, to the model described in |]14| . 



This group recently extended this "screening corrections" model to con- 
sider explicitly gluon emission in the photon wavefunction. For more details 
we urge the reader to consult the recent preprints 

A more precise treatment of the two gluon exchange graphs relates a{x, d±) 
to the unintegrated gluon distribution function /(x, where l± is the trans- 
verse momentum in the loop below the quark box. Since the gluon density 
itself is very poorly constrained experimentally in the small-x region, the 
prospects for accurately constraining f{x,l±), which is related to its deriva- 
tive with respect to scale, from current data alone are very poor. However, if 
one imposes an additional constraint on the unintegrated gluon, for example 
that it obeys the CCFM equation [53 then the distribution becomes more 



constrained (see |Q for a recent attempt to extract f{x,k±) along these 
lines). Unfortunately, this involves making a set of choices on precisely how 
to do this and hence introduces a certain additional model dependence into 
the answer. 

To summarize, each of the four models for a that we discussed contain 
a perturbative piece that rises steeply with energy at small dj_ and a part 
which rises much more slowly with energy that dominates at large transverse 
distances. In HERA kinematics, in each of the models the short distance part 
of a reaches magnitudes typical of those observed in hadronic interactions 
(tens of millibarns). The rapid rise of this piece with energy is tamed in three 



of the models we discussed M, PTSI, IIB 
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Figure 1: The structure function F2{x,Q'^) as a function of Q^. The 1997 
ZEUS data reveals the approximate F2 oc behaviour for the lowest 
values of Q^. 




Figure 2: The logarithmic slope of F2, derived from HI data, plotted as a 
function of x at fixed values of > 3.5 GeV^. A monotonic increase of 
the slope with decreasing x is observed at fixed Q^. The dashed curves are 
extrapolations to regions outside of the data region considered in the fit. 
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Figure 3: The logarithmic slope of F2 derived from ZEUS data plotted as a 
function of x at fixed values of Q^. A monotonic increase of the slope with 
decreasing x is observed at fixed Q"^. 



25 



bl-2 



a> 1 



^0.8 
-o0,6 



0,4 
0,2 




bl.2 
<^ 

o 

^ 1 



-00,6 



0,4 - 
0.2 7 

- 



10 



,-6 



ZEUS Preliminary 



10 



-1 



1^ 



. * « 



• W=210GeV 
W=130GeV 

★ W=85GeV 

* W=50GeV 
■ W=30GeV 
1^ W=14GeV 



1%, * 
I t ■ 



10 



-5 



• 



• 
• 



10 



10 



10 



10 



ZEUS Preliminory 



Q' (GeV') 



• W=21 OGeV 
W=130GeV 

★ W=85GeV 
W=50GeV 
■ W=30GeV 

* W=14GeV 



10 



-4 



10 



-3 



10 



-2 



10 



Figure 4: The logarithmic slope of F2 derived from ZEUS data plotted as a 
function of Q"^ (upper plot) and x (lower plot) for fixed values of W. Note 
the apparent turnover in both plots for fixed W. 
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Figure 5: The profile in d± for tlie pliotoproduction of J/ip, using the dipole 



model of |14|, at three different energies. 
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Figure 6: The ratio of two models for a: the QCD model [11^ divided by the 



saturation model 16 . 
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